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I, INTRODUCTION 
1.1 BACKGROUND. 


In the past several years, highly accurate models of atmospheric 
transport of x-rays and material response of homogeneous and composite 
materials to x-ray energy deposition have been developed at the Ballistic 
Research Laboratory. Calculations have been performed involving x-ray 
transport in both the atmosphere and in targets, energy deposition in 
target materials including ablators, substructure and components, and the 
hydrodynamic material response of structures to energy deposition. 

The effects of various parameters on material response prediction of 
vulnerability and lethality of US ABM's and foreign RV's have been 
determined by relating these effects to critical radii. Most recently, 
models of the nonequilibrium atmospheric response to nuclear bursts have 
been used in parameter studies involving the prediction of radar obscu- 
ration volumes. 


1.2 X-RAY TRANSPORT AND MATERIAL RESPONSE. 


The uncertainties present in the production of the atmospheric trans- 
port data involve the methods of solving the transport equation (e.g. Monte 
Carlo or finite difference), the cross sections used in the calculations, 
and the question of whether atmospheric response can be ignored in the 
calculation of absorption coefficients. These uncertainties produce 
uncertainties in the radii for which various material response phenomena 
begin to occur. 


Calculations have been performed to determine x-ray fluences delivered 
to a target predicted by both a responding medium model and a cold air 


cross section mode1?2, Indications are that while the x-ray free-field 
data obtained using the cold air approximation may not suffice for all 
x-ray/target interactions of interest, particularly at high altitudes, 
the cold air approximation certainly suffices for altitudes below 30 
kilometers and for the co-altitude ranges of interest. 


During FY72 an x-ray ,data base” was generated using the Monte Carlo 
code OGRE, and a code FOX’ was written to fold arbitrary weapon spectra 


1. Joseph Lacetera, "Effects of Nonequilibrium Phenomena Upon X-Ray 
Transport in the Atmosphere (U)'', US Army Ballistic Research 
Laboratories Report No. 1636, SRD, February 1973. AD #524945L 

2. Joseph Lacetera, "Effects of Cold Air Cross Sections on X-Ray Transport 
in the Atmosphere (U)'', US Army Ballistic Research Laboratories (draft). 

oe Norman E. Banks and Wayne A. Coleman, "Transport of Photons Through 
Air Using Source Energy Band Structure from 300 keV to 2 keV (U)", 

US Army Ballistic Research Laboratories Report No. 1577, U, April 1972. 

4. Christopher M. Clayton, "A Users Manual for the BRL Folding Code, 

FOX (U)", BRL Report No. 1655, US Army Ballistic Research Laboratories, 
U, June 1973. AD #912571L 
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into the data base in order to generate real transported x-ray environ- 
ments. The results of FOX folding calculations compared favorably with 
previous calculations using both Monte Carlo and discrete-ordinates 
techniques, and a handbook of x-ray transport data for candidate ABM 


5 
warheads was generated . 


During FY73 the need for angular dependence in the x-ray free- 


‘ ‘ : ri 
field data used in energy deposition calculations was established® : 


Subsequently, during FY74 an energy-, angular-, and time-dependent x-ray 
data base was generated with the newly developed, BRL-sponsored parT® 
code, and the x-ray deposition code FSCATT was medatied” to accept 


angular dependent data. During FY75 this data base was folded with 


weapon spectra newly generated under BRL sponsorship, using FLAIR!® to 


produce angular-~-dependent input for x-ray deposition calculations with 
the new FSCATT. 


During FY74 and FY7S studies were performed to relate uncertainties 
uncovered in previous calculations to their effects on critical radii 


S. Joseph Lacetera and Christopher M. Clayton, ''X-Ray Transport Data 
for Candidate ABM Warheads (U)'', US Army Ballistic Research Labora- 
tories Report No. 1697, SRD-N, January 1974. AD #529019L 

6. Janet Lacetera and Joseph Lacetera, "Sensitivity of X-Ray Deposition 
in ABM Materials to Free-Field Angular Dependence (U)", US Army 
Ballistic Research Laboratories Report No, 1674, SRD-N, September 
1973, AD #527474L 

7. Janet Lacetera and Richard V. Haire, "Sensitivity of X-Ray Energy 
Deposition Profile in ABM Materials to Differences in X-Ray Free- 
Field as Generated by the FOX and ATR Codes (U)", BRL Memorandum 
Report No. 2451, S, March 1975. AD #C001819L 

8. W.A. Woolson, D.L Huffman, W.H. Scott, Jr., W.A. Coleman, 'DART- 
A Monte Carlo Code for Atmospheric Transport of X-Rays and Gamma- 
Rays (U)', BRL Contract Report No. 198, U, December 1974 (AD 
BO01107L. 

9. R.H. Fisher, ''FSCATT: Angular Dependence and Filter Options (U)", 
BRL Contract Report No. 290, March 1976. AD #A023188 

10. W.H. Scott, Jr., and B.L. Colborn, ''FLAIR - A Scaling and Folding 
Code for the Generation of Photon Transport Results in Air", 
Science Applications, Inc., La Jolla, CA, BRL Contract Report 
(to be published). 


(sure-kill and sure-safe distances)", These uncertainties involved 
X-ray transport in the atmosphere and in targets, energy deposition in 
target materials, and hydrodynamic material response of targets to 
energy deposition. Among the parameters evaluated and related to criti- 
cal radii were: (1) differences in x-ray transport models (cold air, 
responding atmosphere), numerical techniques (Monte Carlo, finite 
difference) and inclusion of angular dependence in the x-ray free-field; 
(2) inclusion of the fluorescence process in energy deposition in target 


; 3 . ‘ : 
material! ; and (3) inclusion of angular dependence in energy depo- 
sition calculations. 


A comparison was made of the atmospheric transport of x-rays from 
a unit energy source through air calculated by the folding codes FOX 
and FLAIR, and the radiation transport code THUNDERBALL. (See Appendix A.) 
Observed differences in integral quantities predicted by the codes were 
on the order of 5% or less. These results indicated that the methods on 
solving the transport equation, or choosing the cross section sets, involve 
uncertainties which were small compared to uncertainties involved in the 
material response calculations. It was concluded that, at target distances 
for which the time dependent atmospheric response can be ignored, the 
folding codes produce x-ray free field data which are equivalent to that 
which the best explicit transport codes can produce. (The folding calcu- 
lations take less than a minute on a CDC 7600). However, for smaller miss 
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distances the cold air approximation” used in generating the data bases 
may not suffice for predicting x-ray fluences on BMD targets. 


1.3 NUCLEAR FIREBALL-RADAR. 


The ability to accurately calculate the atmospheric ionization due 
to deposition of x-ray energy from SPRINT interceptors and foreign 
reentry vehicles is critical to the accurate prediction of the nuclear 
environment by radar software systems. Radar systems can be defeated 
by absorption, clutter, multipath, reflection, etc. The ability to 
minimize these effects depends on (1) accurate calculations of the 
atmospheric ionization to create a parameterized data base, (2) obtaining, 
by actual measurements, sufficient data from a particular nuclear burst 


Ii, Joseph Lacetera and Janet Lacetera, "Critical Parameters in X-Ray 
Interactions With ABM/RV Materials (U)', US Army Ballistic Research 
Laboratories, BRL MR 2586, February 1976. AD #CO05358L 

12. Joseph Lacetera, "Uncertainties in the Atmospheric Transport of 
X-Rays Related te Uncertainties in Material Response (U)", Fifth 
Symposium on Nuclear Survivability of Propulsion and Ordnance 
Systems, Menlo Park, CA, October 1975. 

13. Joseph Lacetera and Janet Lacetera, "Sensitivity of X-Ray Depo- 
sition in ABM Materials to the Inclusion of Fluorescence (U)'', US 
Army Ballistic Research Laboratories, BRL R 1861, SRD-N, February 
1976, AD #C005331L 
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(or bursts) in order to (3) correlate the measured data with the data 
base such that the history of the burst can be characterized by computer 
simulation as a function of time. Using this prediction technique 
obsauration volumes for the radar systems can be mapped. The effectiveness 
of the radvar system can be enhanced by avoiding searches within the 
obscuration volumes. However, for such a technique to work, it is 
necessary to have an accurate data base which describes the perturbed 
data base. 


It has been shown? 4 that the atmospheric properties most important 
to such predictions are the free-electron densities and free-electron- 
density gradients in the radar path. (Absorption is an important effect 


for electron densities on the order of 10° particles/cc. The range of 


iv’ te" particles/cc represents the range for which clutter is the 


dominant effect.) These variables are sensitive functions of the 
ionizational structure of the atmospheric gas, both at the point of 
interest and in regions closer to the center of burst from which 
radiation, capable of photodetachment of electrons from negative ions, 
is being emitted. 


. 1 a Penner ie 
Previous calculations with a nonequilibrium, atomic-ion model have 
vn that nonequilibrium species populations have a strong effect on 


ill emission. Subsequent calculations!° with THUNDERBALL have 

dicated that population densities of important absorber/emitters, such 
as NO, can be far from equilibrium at times of interest in the radar- 
blackout problem. Thus the assumption of a fireball in local thermo- 
dynamic equilibrium (LTE) will not always suffice in late-time chemistry 
calculations. In such cases, the detailed ionizational structure of 
the fireball must be calculated by a nuclear fireball code which ex- 
plicitly calculates the nonequilibrium interaction between the radiation 
and the atomic ions and the nonequilibrium chemical kinetics of the 
important molecular absorbers. Such a capability is provided in 
THUNDERBALL. 


Il. THE THUNDERBALL CODE 
2.1 GENERAL DESCRIPTION. 


THUNDERBALL (Thermonuclear Unrestricted Nonequilibrium Detonation 
Evaluator for the Radar Bomblight Problem With All Relevant Processes) 
is a nonequilibrium, radiation-hydrodynamics, photon transport and 
atmospheric response code. Photon and mass transport are calculated in 
spherical geometry and ion population densities are calculated directly 
from the interaction between the radiation field and the responding 


14. Bruce Carney, E.L. Lortie, M.D. Kregel, and F.F. Niles, "Atmospheric 
Deionization Following a 5 MT Burst at 30 km Altitude (U)", 
BRL Report No. 1671, C, September 1973 (AD5278834). 

15. Joseph Lacetera, and Gaelen R. Daum, "Effects of Nonequilibrium 
Phenomena on Radar Transmission-Nonequilibrium Ionization", 
BRL Report (in draft). 


atmosphere. The code also has an option for choosing atomic~ion specie 
densities via the Saha equation, if the user wishes to assume LTE. 
However, molecular specie densities are always calculated explicitly. 


The basic sequence of calculations is as follows. An initial zone 
structure is given as input, along with a description of the weapon source 
and the ambient atmosphere (See Appendix B). The code then determines 
the absorptive and emissive properties of each volume of gas from the 
species populations. A time step is chosen such that stability is 
maintained in the solution of radiation transport, hydrodynamics, species 
densities, and energy partitioning. After the interaction of the radia- 
tion and material field is calculated, the radiation transport and 
atmospheric response are calculated. The macroscopic effects of the 
microscopic processes are used as input to the hydrodynamic calculations, 
which in turn feed back new local conditions. Calculational volumes are 
so chosen that variables, differential in r, vary almost continuously 
from one spatial region to the next, thus ensuring that the calculational 
averages will be truly representative of the condition of the radiation 
and material fields in each volume element. The output of the code 
describes the condition of the atmosphere and the radiation field at 
user-specified instants of time for each spatial cell. 


The code is written in FORTRAN IV, and is maintained in a program 
library using the SCOPE, UPDATE system on a CDC-7600. Memory require-~ 
ment is 160K octal using the 60-bit word size. A restart capability 
is provided by storing all of blank common on a permanent’ file at 
specified time intervals. 


2,2 X-RAY TRANSPORT AND ATMOSPHERIC RESPONSE. 


The problem of determining the detailed effect of nonequilibrium 
phenomena on the atmospheric transport of x rays is extremely complex. 
It involves the interaction of a very intense radiation pulse with a 
rapidly ionizing gas which remains in a nonequilibrium configuration 
for times which are long compared to the time of passage of the pulse. 
This means that ionic species population densities must be calculated 
explicitly without resort to LTE assumptions, and the time steps chosen 
for these calculations must be such that changes in physical quantities 
will be small compared to the value of these physical quantities. 


The important physical processes are: 


Li Photoelectric effect 

eae Auger effect/fluorescence effect 
oe Ionization by electron impact 

4, Bremsstrahlung 


S: Scattering. 


fhe photon transport equation is written 


JE(r,v) 


iT 


+0 (vB) sc S$, (W) - & (2) CE, 


where E(r,v) is the energy flux, S, (nv) is the energy emission flux, 


a(r) is the local absorption coefficient, and c is the speed of light. 
THUNDERBALL calculates photon transport using the finite difference 
method in an Eulerian geometry. A rezoning process is used to provide 
a fine spatial mesh in those regions where physical quantities are in 
the greatest state of flux. 


Spatially, the region of interest is initially divided into 50 
vherical shells of equal differential radius. When the front edge of 
the radiation pulse reaches the outermost shell, the 32 innermost shells 
are rezoned to 16 shells of twice the original differential radius. 
The sixteen additional shells thus obtained are added on to the outer- 
most shells and computations continue. This rezoning process is repeated 
whenever necessary until the problem is completed. 


iiven a distribution of local species populations and local temper- 
ature, the radiation sources per unit volume for a given photon band are 
computed, While these calculations are being made, the various contri- 
butions to the derivatives dn; /dt from radiative recombinations (radiative 


source terms) are being tabulated. The absorption coefficient for each 
photon energy band is then evaluated using an average photon energy. For 
this reason photon bands are always chosen so that the various cut-off 


energies fall at the division points between photon bands. 


When the sources and absorption coefficients have been evaluated 
for the given band, the radiative transport of energy is computed to 
determine where the emitted energy is absorbed and how much is radiated 


way. The species changes dn, due to absorption of photons are next 


evaluated from the computed fluxes and the various partial absorption 
coefficients. Finally the effects of the suprathermal electrons are 
evaluated, This cycle of computations is then repeated for the photon 
band with the next lower photon energies, and so on until all photon 
bands have been covered. The scan proceeds from high to low photon 
energies to ensure that the fluorescence radiation from K-shell 
absorption in lithium-like ions is properly included. Finally, when 
all radiative processes at all photon energies have been evaluated and 
their effects accumulated, the collisional contributions and their 
inverses are added at each position in space. 


2.3 ATOMIC ION MODEL. 


The atomic ion model in the code is based on a photoionization 


. : 17 
cross section set, calculated by Daum! ® using the XSECT code . 
XSECT contains a model of photoionization processes for virtually all 
atomic and ionic species of interest. 


The atomic models currently in use in many codes are so crude that 
their only justification is that they were the best available at the 
time of their development. Positive ions are incorporated into these 


models by extrapolating information about the neutral atom! 8 , or by 


assuming a simple (photon energy) > dependence!’. These models provide 
qualitative solutions, but confidence in the quantitative results must 


be extremely low. As an example the RADCHEM codel? assumed a (photon 


energy)”> dependence for K-shell ionization of an average atom. Ioni- 
zation from the L-shell was ignored as no estimate of its cross section 
was known. Daum calculated cross sections for all the ions of oxygen 
and nitrogen. These cross sections were then folded to yield K- and 
L-shell cross sections for each stage of ionization ot an average atom, 
assumed to be 0.8 nitrogen and 0.2 oxygen. These cress sections were 
found to yield significantly different absorption coefficients (by a 
factor of five) for low energy radiation when compared to the results 
from Richard's model (See Appendix C). At higher energies, the two 
models agreed surprisingly well (within 25%}. Being firmly grounded 

in basic physics, the results with Daum's cross sections may be taken 
as strong evidence that previous calculations using Richard's 
"guesstimates" of the cross sections are satisfactory for energies 

away from the photoionization threshold. The confidence one may place 
on the Richard's-model results is much greater now that confirmation 
has been obtained. Om the other hand, where energies close to the 
photoionization threshold are likely to be important, the XSECT model 
enables one to perform more accurate calculations. 


The XSECT cross sections are weighted appropriately, in THUNDERBALL, 
to an average atom approximation of the atmosphere. In this model the 
atmosphere is assumed to consist of molecules A,, atoms A, and their 


ions aM (M=1-7), along with free electrons e . 


16. Gaelen R. Daum, ''Photoionization Cross Sections for 0, N, A&, and 
Their Ions'', BRL Report No. 1894, 5, June 1976. AD #B012453L 

17. Gaelen R. Daum, "The BRL Photoionization Code XSECT", BRL Report 
No. 1869, U, March 1976. AD #B010267L 

18. W.D. Barfield, G.D. Koontz, and W.F. Huebner, J. Quant. Spectrosc. 
and Radiat. Transfer 12, 1409 (1972). 


i9. P.I, Richards, Summary Report on Investigation of Radiation and 
Chemical Calculations, Tech. Op. Res. Report No. TO-B 62-64 (1962). 
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The photoionization cross sections are then used, employing the 
principle of detailed balance, to compute inverse emissive reaction 
rates. Stimulated emission is included in the detailed balance equa- 
tions for exactness. 


The principles of detailed balance and mass action are used exten- 
sively to ensure that all necessary interrelations between the various 
processes are treated properly. The principle of detailed balance 
requires that in thermodynamic equilibrium the rate of a reaction and 
its inverse are equal. The law of mass action says that the products 
of the concentrations of the species which are present and their rate 
constants are a function of temperature alone. One therefore knows the 
ratio of the concentrations and hence the ratio of the rate constants. 
THUNDERBALL uses these principles to determine inverse cross sections 
which are not otherwise known, and also to establish consistent inverse 
cross section pairs. The rate equations will thus allow thermodynamic 
equilibrium to occur under the proper conditions. 


The code allows local kinetic equilibrium velocity distribution 
among the free electrons. However, excess energy from the photoelectrons 
is accounted for in terms of further ionization and conversion to kinetic 
energy of the average particle. Collisional processes which occur as 
the atomic species seek equilibrium are also accounted for in detail. 


Excited states are ignored on the premise that their effects will 
be small compared to gross changes in species population. 


Collisional Ionization 


alid+ ee wate , ier. 


The ions are assumed to be at rest and in the ground state. Cross 
sections are assumed to be independent of the state of ionization. The 


; ‘ i)+ 
number of ionizations of at) 


of energy E is then, 


per unit volume per unit time by electrons 


M, M. V(E) £(E) @ (E,E’) dE dE’ 


where v(E) is the velocity distribution of the incident electrons of 
energy E, f(E) is their energy distribution, and $(E,E’) the cross 

section for the collision. The energy distribution is assumed Maxwellian, 
and is given by, 


2 oat 4 
M £(E)dE = — Ee , 
. kn 3/2 


The inverse (three body) reaction is proportional to 


2 / 4 / Ta 
M, Mi 4, FCE } BE" ’ ) dk’ de 


and since the reactions are inversely related we can write 


K, = K,/S; (kT), 
ae 
where S; (kT) % Pagers is derivable from the Saha equation. The net 


1 
rate of ionization is then 


2 
RM - EMS Mis 


Collisional Dissociation 


Molecular dissociation processes in the atomic-ion model include: 


Cn Ge Fey (oxygen only) 
A+ Aa SRK 
A, + A, _ A, + 2A 


Free-Free Emission (Bremsstrahlung) 
Bremsstrahlung can occur in the Coulomb field of an ion 


oa Ke se a 


or in the field of a neutral atom 
e +A +e +A+hv 


The heavy particle is assumed to stay at rest, so the photon energy 
comes from the loss in kinetic energy of the electron. 


Free-Free Absorption (Inverse bremsstrahlung) 
hv+ee +A 20 + A° 


Here, there is no restriction on the photon energy, so photons from 
fireball emission are readily absorbed. 
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Stimulated Emission 
mas “hk ae 4K 2 he. 


This reaction, while not contributing significantly to the radiative 
field, must be accounted for in the detailed balance equations. Since 
this emission is coherent with, and in the same direction as, the 
incident radiation, it is treated as negative absorption. 


Photoelectric Absorption and Auger Effect 


i i+]l) - ee 
hv + Gl > i =e hv > I, the ionization energy. 


For K-shell absorption in one- and two-electron ions and L-shell 
absorption in other species, these processes are directly reversible. 
For ions of three or more electrons, the process is not considered 
reversible, The calculations are not corrected for recombinations, and 
there is no corresponding emission term. Besides the highly energetic 
photoelectron, an Auger electron is emitted as the ion decays to fill 
the K-shell hole. The excess energy of these electrons and the L-shell 
photoelectrons is converted to further ionization and kinetic energy 
of the average particle. K-shell absorption in lithium-like ions does 
not produce Auger electrons so fluorescence radiation is calculated, 
in this case, as part of the emission term, 


Scattering 


Scattering of photons with energies above 10-keV by electrons is 
accounted for by using Compton cross sections for free electrons. The 
Klein-Nishina cross section is used to compute the absorption coefficient, 
and scattered photons are treated as emissions. 


2.4 CHEMICAL KINETICS MODEL. 


This model is based on a minimum reaction set for NO production 


given by Sutherland and Zinn’, and is complete in the sense that for 
every reaction included we also include the inverse process. Population 
densities are computed for Ny» 0,, N, O, NO, N,0, 0, » O. , and 0,. The 


D 
model also includes photodetachment of electrons from 0 


2 


attachment to 0,. It is coupled to the atomic-ion model via the reactions: 


and electron 


e +N ##N 
e +9 S:0 , 


20. C. Dexter Sutherland and John Zinn ''Chemistry Computations for 
Irradiated Hot Air", Los Alamos Scientific Laboratory, 
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2.5 ATMOSPHERIC EQUATION OF STATE. 


The basic atmospheric equation of state (EOS) is an ideal gas law 
which uses the average atom approximation to describe the state of ioni- 
zation of the gas. 


2 i)+ 
The molecule Ay» atom A, and its ions a(t) are represented as 


follows: 
A,(Ke1), A(K=2), atl (x=3), ...A*’(Ke9), e° (K=10). 


The gas pressure is calculated as 


where p = b> Nx my (m, = mass of specie K), the fi, are 1/2 the degrees of 
freedom of particles of specie K, and the ERK are the corresponding heats 


of formation. 
2.6 HYDRODYNAMICS. 


Given the particle number densities, N,, of the average atom model 
particles with their atomic masses Ms the conservation laws are 
written as follows: 


Conservation of Momentum 


a(py) 
t 


C aP 
3 Ve V toy) + a> =f 


R 


where p = > Ny ms fp is the net force per unit volume due to the 


absorption and emission of radiation, P is the pressure, and v the flow 
velocity. 


Conservation of Energy 


3E. 


at sie "| (E, tN” (Fv) = -¥ fp -V¢ Qe 


where Qe is the net radiative emission energy flux, so V ° Qe is the net 


power emission per unit volume. 


Conservation of Mass 


oN, 


at 


+ V (Nyy) = Cy 


where Cy represents the creation and destruction of particles of specie 


K due to ionization and recombination, 


Heating due to the collision processes causes hydrodynamic expansion 
and consequent shock wave propagation. 


2.7 CODE INPUT. 


To provide a restart capability all of blank common is automatically 
stored by the code on a permanent I/0 file at specified output times. 
These “runs'' are stored consecutively so that a restart can be made 
from any of these points in time. Handling of this I/O file is as 
described in the following card input description. 


The first two input cards are read in alphanumeric format; all 
other cards are input through the use of NAMELIST statements. 


TITLE: FORMAT (12A6), problem identifier 


REELID: FORMAT (Al0), permanent I/O file identifier written ty the code 
at the beginning of a new I/O file, and checked on restarts. 


NAMELIST START (I/O Information) 
IOB: Number of the I/O file on which blank common variables are stored at 
specified output times to provide a restart capability. On initial 
runs IOB=1 or 2. On restarts, IOB must equal 2. 
IOB1: Number of new 10 file, if desired, on restarts. On initial runs 
IOB1=IOB. On restarts: IOB=1 to write on a new file; IOB=2 
to write on the same file as was read. 
NRUN: The number of runs on the restart IO file. On initial starts 
NRUN=0. On restarts NRUN is specified according to the desired 
restart time, as indicated on previous printed output. 


TIMMAX: The maximum time, in CPU seconds (decimal), that the code will 
use in the calculation. 


NAMELIST PROBLM (Problem parameters) 


IHOB: Altitude in kilometers. 
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ID: Weapon source parameter. ID=O if source is radiating; ID=-1 if 
source is off. 


IE: The number of Eulerian zones used at the current time. Initially IE 
may be any number from 1 to 50. 


IPAR: The number of weapon source spectra used to describe the time 
dependent weapon output. 


RMAX: The distance (cm) that 50 zones will initially cover. The code 
divides RMAX by 50 to get initial zone size. 


TIME: The current problem time. 
YIEL: Weapon source yield (calories) 


NAMELIST WEAPON (Source spectrum) 


WOFLUX (t,E): Weapon output flux (cal/keV) as a function of time 
and photon energy. 


NAMELIST RUN (Control parameters) 

CKTIME: Time for next print edit (sec). 

END: Real time (sec) at which calculation will terminate. 
FREQ: Print edit frequency (sec) 


LMAXT: Number of photon energy bands used in radiation calculations. 


NSS1: = 0 for x ray weapon source; 
1 for parameterized initial conditions. 


NSS2: = 0 for calculations with chemical kinetics; 
= 1 for calculations without chemical kinetics. 
NSS3: = -1 on initial start; = 1 for expanding grid; = 0 if geometry 
is static. 
NSS4: = 0 for nonequilibrium atomic ion model; = 1 for using Saha 
equation; = -1l for cold air calculations. 
NSSS: = 0 for automatic rezone; = 1 for manual rezone. 


2.8 CODE OUTPUT. 


When the code initially receives control from the system it prints 
the computer clock time (T7600) and all of the NAMELIST variables. A 
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page is then printed which is headed by the name of the problem (TITLE), 
the maximum CPU time before cutoff (TIMMAX), and a list of control 
parameters and physical variables that are all defined at the end of 
this section. This is followed by the weapon source and radiation 
field frequency group edit. The remaining printed output is periodic, 
as specified by the FREQ parameter. 


The first page of the periodic output gives information on the 
problem time, current and previous time steps, time step control 
parameters, and the number of zones in the grid. The second page 
lists, as a function of r, mass density, particle density, particle 
velocity, pressure, temperature, internal energy, electron charge per 
cc, and zone dimensions. The fourth page is an edit of the atomic ion 
model, and the fifth page is an edit of the chemical kinetics model. 
Subsequent output lists information on the radiation field, including 
absorption and emission by the gas, as a function of photon energy. 
At the end of each output print (each run) the number of the run and 
the number of the I/O file are printed. 


The following is a partial, alphabetical list of the output variables. 
Others are self-explanatory. 


ABSORPTION: Energy (joules) absorbed from the total radiation field. 
AM+/voL: Number density (uumole/cc) of species in atomic ion model. 


ATMF: Time step control parameter based on reaction rates in atomic 
ion model. 


AVE VEL: Average particle velocity (cm/sec) in zone I. 


CHMF: Time step control parameter based on reaction rates in chemical 
kinetics model. 


DT: Current time step (sec). 
E: Total excess (above ambient) energy (joules) in problem. 


ECHG: Time step control parameter based on rate of deposition of radiative 
energy into gas. 


EKIN: Total kinetic energy (joules) in gas. 


EMISSION: Energy (joules) radiated from zone I. 


ERAD: Total radiation field energy (joules). 


FLCH: Time step control parameter based on radiative flux rates of 
change. 


HYVDR: Time step control parameter based on hydrodynamic calculations. 
MASS DENS: Density (gm/cm?) in zone I. 
NCYC: Number of calculational cycle. 


NUCL DENS: Total particle density (yyumole/cc) in zone I in atomic 
ion model. 


PART/VOL: Total particle density (1/cc) in zone I, in chemical kinetics 
model; particle densities for individual species are fractions 
of this number (e.g. NITROGEN has the value .800 for ambient 
air). 


PRESSURE: Kinetic pressure (dynes /cm”) 


TEMP: Kinetic temperature (K). 


III. SAMPLE CALCULATION 


The configuration of a sample input deck for a THUNDERBALL calcu- 
lation is shown in Figure 1. 


The first card (TITLE) is a 72-character alphanumeric identifier 
describing the problem. The second card specifies the alphanumeric 
identifier (REELID) of the I/O tape (or disk file) on which all of 
COMMON is stored periodically to provide a restart capability. The 
remainder of the input is via the NAMELIST option. All other required 
information is stored in the code. 


The START NAMELIST provides I/O information to the code. NRUN = 0 
means that this is an initial run (not a restart). IOB = IOBL = 2 means 
that there will be one I/O file called TAPE 2. TIMMAX = 60 means that 
the code will stop calculations, with a normal termination after 60 
decimal seconds of CPO time. 


The PROBLM NAMELIST indicates the type of problem being run. 
IHOB = 30 means that the burst altitude is 30 kilometers. ID = 0 normally 
implies a radiating photon source however this control is overridden in 
this particular problem by the RUN variable NSS1. IE = 50 means that 
the problem will initially have 50 spherical shells spanning the 
distance RMAX = SE+0S cm in equal increments. The total yield in the 
calculations will be initiated at TO = 1.0 sec. IPAR is also over- 
ridden in this problem, by the RUN variable NSS1. 
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With the current coding THUNDERBALL expects a description of a 
photon source in all problems. In this problem no late-time radiation 
was included so the value of 0.0 was read into a single element of 
the werpon-output-flux array. 


The RUN variables are the control parameters which specify the way 
in which the problem will be run. NSSl=1 means that parameterized 
initial conditions will be used; NSS2=0 means that chemical kinetics 
will be included; NSS3=-1 must be set this way on initial starts; 
NSS4=0 means that atomic ion populations will be calculated in non- 
equilibrium; NSSS5=1 means no automatic rezoning. CKTIME (sec) is the 
time selected for the first printout (and COMMON dump to TAPE2); and 
FREQ(sec) is the frequency of printing and dumping. END is the time 
at which the calculation will normally terminate if the TIMMAX variable 
has not already terminated the calculation. Normal terminations are always 
accompanied by a printout and tape dump. LMAXT is the number of photon 
energy bands that will be used in transporting photons; these being 
bands 1 through 17. The output corresponding to this input stream is 
listed in Appendix D, along with a glossary of the output variables. 


IV. CONCLUSIONS AND RECOMMENDATIONS 


The current atomic model could be further improved by incorporating 


the results of Missavage, Manson and Daum?! , It is shown there that 


it is not necessary to calculate cross secticns for each stage of 
ionization, as Daum has done, but that it is sufficient to calculate 
cross sections for the isonuclear sequence. This simplification in 

the model represents a great saving of labor and computer time. (Rather 
than doing thirty-four cross section calculations, only twenty-four are 
required.) With these modifications, it is estimated that the atomic 
model portion of the code would run 30% faster. The Manson research 
has proven its applicability and worth to studies of radar obscuration 
volumes, giving greater reliability to the results, allowing further 
simplification of the model, and increasing the range of problems that 
may be treated. The work should prove equally applicable to other 
research of military interest (such as laser propagation), as well as 
outstanding problems in other areas such as plasma physics and astro- 
physics. 
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APPENDIX A: X-RAY TRANSPORT DATA 


A comparison of time-integrated flux (fluence) as a function of 
mass depth has been made between the folding codes FOX and FLAIR and 
the THUNDERBALL code for the same x-ray source. The results are shown 
in Table A as a function of mass depth in arbitrary units. The numbers 
in parenthesis are the ratio of the calculated fluence to that obtained 
using the FOX code. 


TABLE A. 4nr? FLUENCE AS A FUNCTION OF MASS DEPTH 


Mass Depth 
(arbitrary FOX FLAIR THUNDERBALL 
units) 
05 848 (1.00) .922 (1.09) wo2o) (2970) 
.10 .766 (1.00) 7797 (Gle03)) -728 (.940) 
. 30 605 (1.00) ole (1.02) -960 (.950) 
790 .476 (1.00) -491 (1.03) -430 (.900) 
1.00 .268 (1.00) 275 (1.02) ---- 
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APPENDIX B: ATMOSPHERIC DATA 


Atmospheric data used for ambient air were taken from the U.S. 
Standard Atmosphere, 1962. Specifically these data were particle 
density and gas temperature as a function of altitude. Energy density 
in units of erg/yumole was computed from the temperature and installed 
in the code along with particle density in units of uumole/cc. This 
data is shown in Table B, 


| TABLE 8. ATMOSPHERIC DATA 


Altitude Energy Density Particle Density ||| Temperature 
(km) (erg/uumole) (upmole/cc) (°K) 
0 2.533 + 06 4.229 + 07 288 
10 6.624 + 05 1.428 + 07 223 | 
20 1.382 + 05 3.069 + 06 217 
30 2.992 + 04 6.356 + 0S 227 
40 7.179 + 03 1.379 + 05 250 
50 1.995 + 03 3.545 + 04 271 
60 5.614 + 02 1.056 + 04 256 
70 1.380 + 02 3.022 + 03 220 
80 2.594 + O1 6.907 + 02 181 
90 4.110 + 00 1.094 + 02 181 
100 7.537 - Ol 1.727 + Ql 210 
110 1.836 - 01 3.437 + 00 257 
120 6.308 - 02 8.683 - O01 350 
130 3.058 - 02 2.756 - 01 534 
140 1.853 - 02 1.249 - 01 714 
150 1.266 - 02 6.824 - 02 893 


TABLE C.2 COLD AIR ABSORPTION COEFFICIENTS (E, = 1.0 - 300 keV) 


Photon Energy Absorption Coefficients 
(keV) - (om) 
Ref 18 Richards FLAIR 
1.0 - 2.0 2.05 + 03 1.43 + 03 0.0 + 00 
2.0 - 4.0 2.87 + 02 1.79 + 02 0.0 
4.0 - 6.0 4.24 + 01 2.57 + O1 0.0 
6.0 - 8.0 1.34 8.36 + 00 1.04 + 01 
8.0 - 10 5.70 + 00 3.74 6.91 + 00 
iG = ES 2.44 1.65 2.38 
15 - 20 7.30 - 01 5.35 - Ol 8.23 - 01 
20 - 25 3.06 2.40 3.68 
25 = 30 1.55 1.28 1.93 
30 - 40 7.65 - 02 6.68 - 02 8.86 - 02 
40 - S50 Sab? 3.00 ; 3.91 
50 - 60 1.60 160 2.03 
60 - 70 9.06 - 03 9.58 - 03 at? 
70 - 80 5.60 6.18 7.36 - 03 
80 - 90 3.69 4.22 4.88 
90 - 105 2.39 2.84 3.07 
105 = 120 1.48 1,83 1.93 
120 - 140 9.36 - 04 1.20 1.22 
140 - 160 5.82 7.73 - 04 7.66 - 04 
160 - 190 3.60 4.95 4.59 
190 ~- 220 2.13 3.04 2.78 
220 =- 260 1.30 1.91 Lvs 
260 - 300 7.84 - 05 1.19 1.07 
30 


APPENDIX D: SAMPLE OUTPUT 
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